Spin squeezing: 
transforming one-axis-twisting into two-axis-twisting 
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Squeezed spin states possess unique quantum correlation or entanglement that are of significant 
promises for advancing quantum information processing and quantum metrology. In recent back to 
back publications [C. Gross et al, Nature 464, 1165 (2010) and Max F. Riedel et al, Nature 464, 
1170 (2010)], reduced spin fluctuations are observed leading to spin squeezing at — 8.2dB and — 2.5dB 
respectively in two-component atomic condensates exhibiting one-axis- twisting interactions (OAT). 
The noise reduction limit for the OAT interaction scales as cx 1/N 2 ^ 3 , which for a condensate with 
N ~ 10 3 atoms, is about 100 times below standard quantum limit. We present a scheme using 
repeated Rabi pulses capable of transforming the OAT spin squeezing into the two-axis-twisting 
type, leading to Heisenberg limited noise reduction oc 1/N, or an extra 10-fold improvement for 
N ~ 10 3 . 

PACS numbers: 42.50.-p, 03.75.Gg 



Squeezed spin states (SSS) 0, 0] are entangled quan- 
tum states of a collection of spins in which the correla- 
tions among individual spins reduce quantum uncertainty 
of a particular spin component below the classical limit 
for uncorrelated particles Q . Research in SSS is a topical 
area due to its significant applications in high-precision 
measurements [3l4l0l| and in quantum information science 
Ill4l5|. Squeezed spin states were first introduced by 
Kitagawa and Ueda, who considered two ways to pro- 
duce them. The simplest to implement uses a "one-axis 
twisting" (OAT) Hamiltonian, but the state it produces 
does not have ideal squeezing properties. A more com- 
plex approach uses a "two-axis twisting" (TAT) Hamil- 
tonian and produces an improved state. Other mecha- 
nisms for producing SSS have also been investigated, es- 
pecially those based on atom-photon interactions r| 12 1 
and quantum non-demolition measurements 



Atomic Bose-Einstein condensates are promising sys- 
tems for observing spin squeezing. Assuming fixed spa- 
tial modes, condensed atoms are described by a collection 
of pseudo-spin 1/2 atoms, with spin up (|t)) and down 
(|j)) denoting the two internal states or spatial modes 
12,22- 25| . The two recent experiments [J Q raise sig- 
nificant hope for reaching the theoretical limit of spin 
squeezing cx 1/N 2 ^ 3 with N the total number of atoms 
for the OAT model Q. Both experiments utilize two in- 
ternal hyperfme states of condensed atoms, with the OAT 
interaction cleverly constructed from binary atomic col- 
lisions, possibly accompanied by systematic and funda- 
mental imperfections not confined to the two state/mode 
approximation. They can be further degraded by atomic 
decoherence and dissipation 0, [f| . 

This Letter describes a readily implcmentable idea for 
improved spin squeezing in the two experiments. Given 
the reported OAT model parameters 0, Q , we propose 



a coherent control scheme capable of transforming the 
OAT into the effective TAT spin squeezing, leading to a 
Heisenberg limited noise reduction cx 1/N, or a further 
10 fold improvement for a condensate with ~ 10 3 atoms. 

With a coherent coupling, the two-component conden- 
sate is described by the following Hamiltonian 



H = X J 2 z +mJy, 



(1) 



as in 0, Q . The collective spin J is defined according 
to J v = J^fcOv /2 in terms of the Pauli matrices a*, 
(v = x, y, z) for the pseudo-spin of the fc-th atom. The 
first term on the right-hand side of Eq. ([TJ is the nonlin- 
ear interaction responsible for the OAT-type spin squeez- 
ing, with x the atomic interaction parameter. The Rabi 
frequency f2(f) = Qof(t) results from near- resonant two- 
photon (microwave+rf) coupling between the hyperfme 
states 0, [H[ ■ Rs maximum amplitude £1$ and the tempo- 
ral envelop f(t) can be controlled experimentally 0, 

Our idea for transforming the OAT {TJ into the TAT 
makes use of multiple 7r/2 pulses affected with the cou- 
pling term HJ y . In the Rabi limit, |f2| 3> iV|x|, nonlinear 
interaction can be neglected while the collective spin un- 
dergoes driven Rabi oscillation. A pulse with an area of 
7r/2 corresponds to J_ Cl(t)dt = n/2, which gives the 
transformation 
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tt/2 



(2) 



where Rg = e~ », rotating the spin around the y-axis 
by an angle 9. The proposed pulse sequence is periodic 
as shown in Fig.QJa). Each period [Fig.QJb)] is made up 
of the following: a 7r/2 rotation about j/-axis (red pulse), 
a free evolution for 2St, a second — 7r/2 rotation about 
y-axis (blue pulse), and a second free evolution for St. 
The period is t c ps 35t, neglecting the time needed for 
affecting the two ±7r/2 pulses. 
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FIG. 1: An illustration of the proposed pulse sequence with 
pulse envelop Q/Qo vs. time t in arbitrary units, and the 
associated spin distributions on Bloch spheres, (a) The multi- 
pulse sequence is cyclic, with one period between nt c and 
(n + l)t c (shaded) and zoomed in (b). Each period contains 
two 7r/2 pulses. A red pulse rotates the spin from along the 
positive z-axis to the ir-axis, and a blue one coherently phased 
to achieve the opposite, is a —tt/2 pulse, as shown on the large 
Bloch sphere of (b). The four small Bloch spheres are located 
at their corresponding times. The spin distributions for the 
upper two are centered on the z-axis, while the lower two are 
centered on the z-axis, due to the applications of red and blue 
7r/2 pulses. The twistings from the OAT interaction (fl} are 
illustrated with counter-rotating circular arrows pierced by 
the x- or z-axis of all Bloch spheres. 



The time evolution operator at time t = nt c (n - 
1, 2, 3, ■ • • ), i.e., after n periods, is given by U n , where 



U 



t2 



Z p „ — 2irJ p — irJ z — 2irJ„ 



(3) 



is for a single period and r = Using Baker- 

Campbell-Hausdorff formula, we find 

U = e-^ 2J *+ J * 2 ) exp[*T 2 {J x , {J v , J z }} + 0(t% 

where {A, B} = AB + BA denotes anti-commutator for 
operators A and B. Expanding for small r, we arrive at 
exp[iT 2 { J x , { J y , J z }}+0(t 2 )] ~ 1 after neglecting higher- 
order terms. Hence we end up with U ~ e 
and the time evolution operator U n : 



-inr{2Jl+Jl) 



-i X {2Jl+Jl)t/S 



equivalent to that given by an effective 



Hamiltonian H eS = x(2J 2 + J 2 2 )/3 = x(J x ~ Jy 



J 2 )/3. 

After dropping a constant J 2 = j(j + 1), H e g reduces to 
that of the TAT. This shows the application of our pro- 
posed 7r/2 pulse sequence effectively transforms the OAT 
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of the recent experiments into the TAT mod- 
prtjividcd r <C (2N)~ [26[. For a fixed Xi t ne time 



TAT 



the 



ptimal spin squeezing from the TAT cx x( J 2 — J 2 ) is 
L/3 of that from the OAT cx %J 2 with N = 1250, 
squeezing occurs around 3 times faster for the TAT. 

iespite of the three times reduction in the effec- 
stiength (x/3) of the transformed Hamiltonian H e g, 
tiixe for observing the optimal SS remains almost the 
Consequently, we expect degradation from atomic 
will be similar to the case of OAT [27j . 
validity for our idea of transforming the OAT into 
can be directly checked through comparing the 
e dynamics from H c h with the actual dynamics 
{T]) accompanied with the sequence of pulses. For 
purpose, we expand the state vector time evolved 
initial state \j,j), 



- iM - J y )/3 \j,3) 



Cm(t)\j,m), (4) 



into eigenstate of J z with j = N/2. The primed summa- 
tion implies m = — j, — j + 2, • • • , j, or m = —j + 1, — j + 
3, • • • , j, for even or odd N respectively 0]. 

First we consider the case of small TV, e.g., TV = 2. 
The quantum dynamics of the effective TAT Hamilto- 
nian H e ff = x{Jx ~ Jy)/3 give rise to time evolved prob- 
ability amplitudes c_i(i) = — isin(%t/3), co(t) = 0, and 
c\(t) = cos(xT/3). For each period t c = 3St of the ex- 
act dynamics given by Eq. ([3]), the probability ampli- 
tudes become c_i(t) = e -^( t +» t <=)/2, c {t) = e~ ix5t /V2, 
Cl (t) = e -ix{t+nU) j2 f or ( n _ i) tc < t < ( n -l)t c + 25t, 
and c_i(t) = -ie^ l x(*-"* c / 3 )sin(xnt c /3), c (t) = 0, 
d(t) = e -^{t-nt c iz) cos ( xn t c /3) f or nt c - St < t < 

nt c . At t = nt Cl i.e., after n periods, they become 
c_i(t) = -ie- axt / 3 sm( X t/3), c (t) = 0, and a{t) = 
e -«2xt/3 cos ( x ^/3) ; exactly the same as that from the ef- 
fective TAT dynamics, apart from an overall phase for 
the J 2 term. 

Next we consider larger TV numerically, e.g., TV = 1250 
atoms and for different number of pulses TV C during the 
time of optimal squeezing x t ~ 31n(4TV)/(2TV) The 
TAT Hamiltonian H e g can generate optimal spin squeez- 
ing when applied to the initial state = It) 8 0]. 
The state expansion Eq. (U) involves basis \j, m) with 
even or odd m because of the conserved parity exp(z7r J z ). 
As a result, the mean spin is always along the z-axis, i.e., 
(J) = (0, 0, (J z ))- For any spin component normal to 
(J), J-y = J x cos"f + Jj,sin7 with an arbitrary angle 7, 
its variance is found to be 

(AJ 7 ) 2 = -[C + -4cos(2 7 ) + B 011(27)], (5) 



where A = (J 2 - J 2 ) = Sx^M^ff), B = (J x J y + 
J V J X ) =Im(J 2 ), andC = (J x + Jy) = j(j + l)-(J z ). The 
optimal squeezing angle 7 op is obtained from minimizing 
(AJ 7 ) 2 with respect to 7, yielding tan(27 op ) = B/A. As 
the coefficient A — Rc(J'^_) = 0, we have 7 op = 7r/4 or 
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FIG. 2: Spin squeezing parameter at different pulse numbers 
N c — 10 (upper), 100 (middle), and 1000 (lower) compared 
with the effective TAT H c ff dynamics in black solid lines with 
the actual ones (0 using the proposed sequence of pulses in 
red solid lines. The black dot-dashed lines denote OAT re- 
sults, all for N — 1250 atoms. 



37r/4, which implies that squeezing and anti-squeezing 



occurs for J x ± J y , 



along the angle bisector of x- 



and y- axis, with the reduced variance V- = (C — \B\)/2. 
Alternatively, the increased variance V+ = (C + \B\)/2 is 
associated with anti-squeezing. 

The reported degree of squeezing is measured by spin 
squeezing parameter Q = NV-/ (J) 2 from Ramsey spec- 
troscopy [3j, which differs slightly from £ 2 = 2V-/j. The 
latter form is used to graph spin fluctuations in the rel- 
evant figures of Rcfs. [j, which is independent of 
angular momentum coordinate system or specific mea- 
surement scheme. For the coherent spin state the 
variances (AJ 7 ) 2 = j/2 and £ 2 — 1. Time evolution 
from H c ff, develops quantum correlation, which trans- 
forms \j,j) into a SSS with £ 2 < 1. As shown in Fig. [5J 
with the increase of pulse number N c , the actual OAT 
dynamics approach and eventually settle down to the ef- 
fective dynamics of H c r . The minimum number of peri- 
ods required before reaching the optimal squeezing point 
is found analytically and checked numerically to satisfy 
N c ^> ln(4AT) [26| in order to generate the effective dy- 
namics of the TAT. For N = 1250, ln(4JV) « 8.5, thus 
it is reasonable to expect a near perfect agreement at 
N c = 1000. In practice, without special control tech- 
niques, errors from repeated pulses can build up, limiting 
the performance of our proposal. 

Fortunately, however, we find that even for signifi- 
cantly fewer number of pulses, the TAT-type spin squeez- 
ing limit can still be reached at selected times before ar- 
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FIG. 3: Top-higher panel: Spin squeezing parameter calcu- 
lated from the actual dynamics of (JXJ) for N c = 20 pulses 
shown by red and blue dots respectively for time windows af- 
ter the first (red) and second (blue) pulses of each cycle. For 
comparison, the black solid (dot-dashed) line denotes results 
from the associated TAT (OAT) model. The top-lower panel 
shows the corresponding optimal squeezing measurement an- 
gle. The dichotomy of two narrow distributions around ±7r/4 
is due to the repeated n/2 pulses rotating between the x- and 
z-axis. For the TAT around the z- or i-axis, this angle is 
simply a flat distribution at ±tt/4, as shown by black dashed 
lines. At the reduced number of cycles N c = 20, the TAT is 
not yet fully effective, which is the cause for small slopped 
distributions through ±7r/4. The slops eventually reduce to 
zero with the increase of iV c . The bottom panels display in 
detail the shaded region where the slops are now displayed in 
absolute values, all for TV = 1250 atoms. 



riving at the optimal squeezing. This is presented in the 
upper panel of Fig. |3l where spin squeezing obtained 
from the actual dynamics ([1]) are shown by red and blue 
dots respectively for times after the first (red) and second 
(blue) pulses of each period. One can see clearly that in 
segments of red and blue dots, the oscillating squeezing 
parameter kisses the TAT spin squeezing represented by 
the black solid curve repeatedly. Furthermore, from a 
technical point of view, the measurement of spin squeez- 
ing can be carried out more straightforwardly as the angle 
for optimal squeezing 7 op is essentially fixed, corresponds 
to a nearly flat distribution in the immediate temporal 
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neighborhood as shown in the lower panel of Fig. [3l 

The dynamic behavior of spin squeezing from the 
Hamiltonian (TTJ) has been studied before [24I I25I 1 , in clud- 
ing the idea of turning off the OAT interaction [28[ , and 
the projection into a desired form [29[ as in NMR dy- 
namic decoupling. Our idea of multiple pulses as pre- 
sented in this work is straightforward. It comes with a 
clear physical picture, and makes use of coherent con- 
trol techniques [5(1 E|. With suitable generalization, 
our work suggests that for systems of identical spin-1/2 
atoms, any forms of binary interactions, e.g., the famous 
Lipkin-Meshkov-Glick model 32(, can be transformed 



into the effective TAT-type provided their original spin- 
spin interaction is NOT SU(2) symmetric, or NOT pro- 
portional to J 2 = J 2 + Jy + J%. Wc can further adopt 
techniques like spin echo using additional 7r-pulscs middle 
way between all neighboring ±7r/2-pulses, to effectively 
suppress both spin-spin relaxation or other inhomoge- 
neous effects. 

Before conclusion, we note that in all our simulations 
for the actual dynamics of ([T]), we have assumed the 
±7r/2 control pulses are executed quickly, during which 
the nonlinear twisting interaction is negligible. The re- 
sults presented above in Figs. [2] and [3] are found to re- 
main unchanged irrespective of whether square shaped or 
Gaussian shaped pulses are used as long as they give the 
correct pulse area of it/ 2 and are short enough so that 
spin squeezing during the pulses can be neglected. When 
the phase accumulated from spin squeezing interactions 
are included, they generally lead to a gradual degrada- 
tion of the achievable spin squeezing. Based on numerical 
simulations over a broad range of parameters 0, HH , we 
conclude for large N (> 100) the time for optimal squeez- 
ing yt nn ^ 1.581n(iV)/iV [26j, consistent with the earlier 
result [34|. 

In conclusion, we present an idea based on coherent 
control theory. Adopting repeated Rabi pulses we trans- 
form the OAT spin squeezing observed in the two recent 
experiments 0, @ into stronger and more effective TAT 
spin squeezing. Both analytical analysis and numerical 
simulations are presented that confirm the validity of our 
proposal. The Hcisenberg limited noise reduction oc 1/N 
from our proposal is an additional factor oc I/Y 1 / 3 lower 
than that reached by the OAT model Q. This work 
thus enables improved spin squeezing from oc I/Y 2 / 3 to 
oc 1/N despite a factor of 1/3 reduction in the effective 
atomic nonlinear interaction strength. Additionally, the 
optimal squeezed direction/angle for the TAT interaction 
is free from the swirling associated with the OAT [1, 0, [l| , 
allowing for simpler and cleaner measurements. Our key 
idea of controlled dynamics to affect a TAT interaction 
Hamiltonian from a OAT form can be further general- 
ized to other s yste ms where OAT spin squeezing are dis- 
cussed [l2, 35j 36 1 . Our proposal calls for no major com- 
plications to the available experimental setups, except 
for a more advanced timing sequence for executing Rabi 



pulses while keep tracking of their phases. Thus we fully 
expect that our idea can be realized within current ex- 
periments and hopefully implemented immediately The 
experimental demonstration of our proposal could sig- 
nificantly push the frontier of quantum metrology with 
squeezed spin states into new territory. 
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